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1.  Introduction: 


During  the  contract  period,  the  principal  investigator  H.J.S.  Fernando,  graduate  students 
Rajka  Krstic  and  Eric  Strang,  post-doctoral  fellow  Andrew  Folkard,  and  visiting  scientists 
Professors  Eliezer  Kit  (Tel  Aviv  University)  and  J.C.R.  Hunt  (U.K.  Meteorological  Office) 
worked  on  several  problems  related  to  turbulent  mixing  in  stably  stratified  fluids.  In  addition,  they 
initiated  an  experimental  program  on  the  coastal-ocean  bottom  boundary  layer.  The  results  of  the 
completed  studies  are  outlined  below.  The  experiments  on  the  bottom  boundary  layer  are  still  in 
progress  as  a  part  of  a  continuation  grant.  Relevant  publications  originated  from  the  P.I.'s  group 
during  the  contract  period  are  listed  at  the  end  of  this  report. 

Measurement  of  Turbulence  Near  Shear-Free 
Density  Interfaces 

An  experimental  study  was  carried  out  to  investigate  the  structure  of  turbulence  near  a 
shear-free  density  interface.  The  experimental  configuration  consisted  of  a  two-layer  fluid  medium 
in  which  the  lower  layer  was  maintained  in  a  turbulent  state  by  an  oscillating  grid.  The 
measurements  included  the  root-mean-square  (r.m.s.)  turbulent  velocities,  wave-number  spectra, 
dissipation  of  turbulent  kinetic  energy  and  integral  lengthscales.  It  was  found  that  the  introduction 
of  a  density  interface  to  a  turbulent  flow  can  strongly  distort  the  structure  of  turbulence  near  the 
interface  wherein  the  horizontal  velocity  components  are  amplified  and  the  vertical  component  is 
damped.  The  modification  of  r.m.s  velocities  is  essentially  limited  to  distances  smaller  than  about 
an  integral  lengthscale.  Inspection  of  spectra  shows  that  these  distortions  are  felt  only  at  small 
wave  numbers  of  the  order  of  integral  scale  and  a  range  of  low-wave  numbers  of  the  inertial  sub¬ 
range;  the  distortions  are  pronounced  as  the  interface  is  approached.  Comparison  of  the  horizontal 
velocity  data  with  the  Rapid  Distortion  Theory  (RDT)  analyses  of  Hunt  and  Graham  (J.  Fluid 
Mech.,  84,  209-235,  1978)  and  Hunt  (J.  Fluid  Mech.,  138,  161-184,  1984)  showed  a  qualitative 
agreement  near  the  interface  and  a  quantitative  agreement  away  from  the  interface.  The  RDT 
predictions  for  the  vertical  component,  however,  were  in  general  agreement  with  the  data.  The 
near-interface  horizontal  velocity  data,  however,  showed  a  quantitative  agreement  with  a  model 
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proposed  by  Hunt  (1984)  based  on  non-linear  vortex  dynamics  near  the  interface.  The  effects  due 
to  interfacial  waves  appear  to  be  important  for  distances  less  than  about  10%  of  the  integral 
lengthscale.  As  a  result  of  the  non-zero  energy  flux  divergence,  the  introduction  of  a  density 
interface  to  oscillating  grid  turbulence  increases  the  rate  of  dissipation  in  the  turbulent  layer  except 
near  the  interface  where  a  sharp  drop  occurs.  These  measurements  provide  useful  information  on 
the  structure  of  turbulence  in  shear-free  boundary  layers,  such  as  atmospheric  and  oceanic 
convective  boundary  layers,  thus  improving  modeling  capabilities  of  such  flows. 

Frequency  Spectra  in  Zero-Mean-Shear  Turbulence 

Eulerian  frequency  spectra  were  measured  in  zero-mean-shear  (oscillating-grid  induced) 
turbulent  flows  and  were  compared  with  the  spectral  form  proposed  by  Tennekes  (J.  Fluid  Mech., 
67,  561,  1975).  The  aim  was  to  verify  some  of  the  available  models  for  the  four-dimensional 
space-time  spectrum  of  homogeneous  turbulence.  This  specturm  is  of  immense  utility  for  the 
Rapid  Distortion  Theory  calculations  of  turbulence  near  density  interfaces.  A  good  agreement 
between  the  theoretical  prediction  and  experimental  results  were  obtained  in  a  limited  frequency 
range.  The  empirical  constants  pertinent  to  the  spectral  law,  obtained  experimentally,  were  in  good 
agreement  with  numerical  simulation  results. 

Evolution  of  Kelvin-Helmholtz  Billows 

A  mixing  mechanism  prevalent  in  natural  flows  is  the  formation  and  breakdown  of  vortical 
billows  known  as  Kelvin-Helmholtz  (K-H)  instabilities.  Laboratory  experiments  were  carried  out 
to  study  certain  key  features  of  K-H  billows,  wherein  the  billows  were  generated  in  a  two-layer 
stratified  tilt-tank.  It  was  shown  that  small-scale  turbulent  mixing  is  present  within  billows  from 
the  early  stages  of  their  evolution,  but  mixing  becomes  intense  and  the  billows  are  destroyed  as 
they  achieve  a  maximum  height  and  initiate  collapse  at  a  non-dimensional  time  of  AUt  /X  ~  5, 
where  AU  is  the  velocity  shear  and  X  is  the  wavelength.  When  AUt  /  X  <  5,  the  Thorpe  scale  Lt 
and  the  maximum  Thorpe  displacement  (LT)max)  normalized  by  the  local  billow  height  Lb,  are 
independent  of  both  the  horizontal  location  within  the  billow  and  time  with  Lx/Lb  « (0.49  +  0.03) 


3 


and  (LT)max/Lb  « (0.89  +  0.02).  After  the  collapse  starts,  however,  the  pertinent  lengthscale  ratios 
in  the  “core”  of  the  billow  show  values  similar  to  those  of  fully-developed  turbulent  patches,  i.e., 
Lx/Lb  ~  (0.29  ±  0.04)  and  (LT)max/Lb  ~  (0-68  ±  0.04).  The  field  observations  were  found  to  be 
in  good  agreement  with  laboratory-based  predictions. 

Resuspension  of  a  Particle  Bed  by  a  Round  Vertical  Jet 

An  experimental  study  was  carried  out  to  investigate  the  resuspension  (fluidization)  of  a 
loosely  held  monodispersed  particle  bed  by  a  momentum  jet  discharging  from  below.  The  work 
was  motivated  by  its  applications  to  limnological  situations  where  the  bottom  sediments  are 
suspended  by  a  series  of  jets  fed  by  a  groundwater  karstic  system.  Two  different  flow  regimes 
were  identified,  and  the  conditions  under  which  they  occur  were  mapped  on  a  regime  diagram 
between  the  two  important  non-dimensional  variables,  ho/dp  and  ho/(Mo^^/ws),  where  ho  is  the 
height  of  the  particle  layer,  Wg  is  the  particle-settling  velocity  and  Mo  is  the  momentum  flux  of  the 
jet.  The  maximum  height  of  rise  of  particles  within  the  jet  was  also  determined  as  a  function  of 
these  non-dimensional  variables.  Comparisons  of  laboratory  and  field  observational  results  were 
also  made. 


Turbulent  Wakes  of  Stratified  Flow  Past  a  Cylinder 

Laboratory  measurements  were  carried  out  to  investigate  the  evolution  of  a  turbulent  wake 
behind  a  right  circular  cylinder  moving  in  a  linearly  stratified  fluid.  The  flow  field  is  determined  by 
the  internal  Froude  number  Fi  and  the  Reynolds  number  Re,  but  at  high  Re,  Fi  becomes  the  only 
governing  parameter.  Measurements  show  that  stratified  turbulent  wakes  can  be  classified  into 
three  flow  regimes,  based  on  Fi.  When  Fi  2,  the  wakes  do  not  grow  downstream,  and  remain  at 
approximately  constant  height.  For  2  ^  Fi  3,  the  wakes  grow  to  a  maximum  height  at  M  =5 
and  then  collapse  physically;  for  M  ^  3,  the  maximum  height  is  achieved  at  Nt  ~  2.5,  before  the 
collapse  begins.  The  evolution  of  such  other  length  scales  as  the  Ozmidov,  Kolmogorov, 
overturning  and  Thorpe  scales  and  the  maximum  Thorpe  displacements  were  measured,  and  their 
behavior  in  the  above  Fi  ranges  delineated.  Length  scale  diagrams  for  the  evolution  of  stratified 
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turbulence  in  cylinder  wakes  were  constructed,  and  were  compared  with  previous  theoretical 
predictions.  These  results  provided  new  insights  into  the  evolution,  collapse  and  two- 
dimensionalization  of  stratified  turbulent  flows. 


Migration  of  Density  Interfaces  Subjected  to 
Differential  Turbulent  Forcing 

A  laboratory  study  was  performed  on  the  migration  of  sharp  density  interfaces  in  the 
presence  of  differential  turbulence  levels  across  them.  It  was  shown  that,  in  a  certain  parameter 
range,  the  net  migration  rate  of  the  interface  can  be  parameterized  (to  the  first  order)  as  the  resultant 
of  the  entrainment  rates  based  on  either  one  or  the  other  layer  is  non-turbulent  (i.e.,  stirring  is  from 
one  side  only).  The  interface  initially  migrates  in  a  preferential  direction  determined  by  the  r.m.s. 
velocities  and  the  integral  lengthscales  of  turbulence  near  the  interface,  and  then  achieves  a  quasi¬ 
stationary  state  wherein  the  system  runs  down  to  a  homogeneous  state  due  to  the  buoyancy 
transport  through  the  interface. 
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